
QR<-matrix(runif(p*p),nrow=p)
Q<-qr.Q(qr(QR))
D<-diag(rexp(p,rate=1))
A<-Q%*%D%*%t(Q)

